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Abstract 

I ^ ■ A stability analysis is made for a non-singular pre-big-bang like cosmological model 

based on 1-loop corrected string effective action. Its homogeneous and isotropic so- 
lution realizes non-singular transition from de Sitter universe to Friedmann-like uni- 
verse, via super inflation phase. We are interested in whether the non-singular nature 
I ' of the solution would be stable or not in more general inhomogeneous case. Per- 

turbative analysis is made for scalar, vector, and tensor linear perturbations, and 
\jr^ • instability is found for tensor- type perturbation. 

\o : 

1 Introduction 

' One of the most remarkable predictions of General Relativity is the existence of singularities. Particularly 

^\ , at the center of the black hole and in the very early universe, one might interpret the singularities as a 

""^ ' result of the breakdown of General Relativity. If so, the behavior of the system "near" the singularity 

should be explained by some more general theory, preferably including GR as an effective theory in 
particular situation. In this sense, searches for singularity-free cosmological models have been made 
(3JT), by many authors ||. Recently non-singular cosmological models based on superstring theory Q were 

presented, and among these, we will here investigate the model presented in pjfl , which is derived from 
• i-h , the superstring effective action with 1-loop string correction pT| 

XI ■ 



S = J d^x^—gi-R - ~(D$y - -(DaY + -[A x e* - A 2 £(a)]i^ B } (1) 

where R ,<f> and a are the Ricci scalar curvature, the dilaton, and the modulus field, respectively. The 
Gauss-Bonnet curvature is 

R 2 gb = R^ kX R^,x - iRTR^ + i? 2 , (2) 
and is a function of modulus field expressed with Dedekind rj function: 

£(a) = -l n [2eV(*e CT )]. (3) 

Since the non-singular nature essentially depends on the behavior of modulus field a, we will concentrate 
on the modulus field only and actually use the effective action appearing in fl3| : 

S = J dV=5{iiJ - \{Dyf ±Z(<p)Rgb}- (4) 

From this action we will derive the back ground equations of motion and see the behavior of the solution 
in the next section. In the section 3, the method of perturbation is explained and equations of motion 
for perturbative variables are derived, and also the stability of this cosmological model is discussed. The 
results are summarized in the section 4. 
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2 Model and the background solution 

We will assume the homogeneous and isotropic metric with conformal time 77, 

ds 2 = a^fi-drj 2 +^ ij dx i dx j ), (5) 

Ttj = z ; 1 y, 9 ; 5~~ oT^y (^) 
1 + j!C(x z + y 2 + z 2 ) 

where /C = 0, +1, —1 for flat, closed, open universe, respectively. From the action (|4|) one can derive the 
equations of motion as 

y/ 2 =6(H 2 +/C)(l- At#) (7) 

(«' + + /C)(l - + (7Y 2 + /C)(l - = (8) 

<p" + 2H<p' + ^(H 2 + K)H'U = (9) 



2a 2 



Here, the conformal Hubble parameter 



a 



H:=— (10) 



a 



is used and prime (') denotes differentiation with respect to conformal time 77. We denote the derivative 
with respect to the physical time t as dot ('), and usual Hubble parameter as H := -. Since the action 
(U) is invariant under the change of ip's sign, we can choose the initial condition of ip so that p increases 
in time (at least near the initial point). Using 

i = \y (ii) 

instead of (^), which is a good approximation near ip = 0, the solutions of the equations of motion which 
continue to Friedmann like phase (H > 0, H < 0) at t — are shown in the fig.l. In these solutions ip 
increases monotonically, so larger ip corresponds to later time. As can be seen from the figure, there are 
two classes of such solutions — singular and non-singular solutions. Solutions a and b in the fig.l are 
singular at <p = 0, whereas c, d, e and f continues beyond the initial singularity. One may also notice that 
these non-singular solutions approaches de Sitter like solution (H ~ const ^ 0) as ip — > —00, t — > —00. 
This is peculiar to the form of the potential (^Tj) . If we use (^) , the behavior of the solutions near tp = 
are quite similar to those shown in the fig.l, but H finally approaches zero in the infinite past. Although 
the asymptotic behavior is quite different, we will use ( [TT| ) as £ in the actual numerical calculation since 
it is simple, and we are only interested in the behavior of the solution near the singularity. 



3 Perturbative analysis 
3.1 Method of perturbation 

We will consider perturbation of our model(|j) to analyze the stability of the non-singular nature. To 
avoid the gauge ambiguity, we here use the gauge-invariant perturbation method ]T^| . The metric is 
decomposed into the back ground part and the perturbed part, and the perturbed part into scalar-, 
vector-, and tensor-perturbation part: 

9nv = (0) £W + 5 s + S v g^ + S T g^. (12) 

Also, the scalar field (p (or the potential £) is decomposed into its background and perturbed part: 

p = Mip + Sip, 

e = (0) e + ^. (13) 
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Figure 1: Behavior of Hubble parameter versus scalar field <p. Time flows from left to right since in 
these solutions <p is always positive. There are two classes of solutions: singular solutions (a and b), and 
non-singular solutions (c,d,e,f). We can see that the initial singularity is avoided when ip crosses zero. 
For detailed discussions sec |Ol. 



Change of perturbed values under gauge transformations is given by Lie derivatives as: 

6Q^SQ + £ X ^Q. 



(14) 



We will introduce such combinations of perturbed variables that the effects of gauge transformation 
cancel. These gauge invariant variables are used to describe the dynamics of the perturbation. Equations 
of motion for perturbed variables are obtained as perturbed Einstein-like equation: 



(15) 



Each term is also decomposed into scalar, vector, and tensor part. Perturbed Klein-Gordon equation is 
not necessary since it is derived from (15). 



3.2 Scalar part 

Scalar part of the perturbed metric is written in terms of 4 scalar functions (f>, ip, B, and E, 



S S g f *» = a 2 (ri) 1 



B 



(16) 



and following combination of gauge-dependent variables are shown to be unchanged under gauge trans- 
formations: 



$ = 



-\a{B-E' 



q, = </, (B-E'), 

a 

6ip M = 6<p + W(f/(B-E'), 
5^ = <5£ + (0) £'(£ -E'). 



(17) 

(18) 

(19) 
(20) 
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Introducing a variable 



: = * - — 7 f -( H- + /C)^ (9i) - H?*}, 
scalar part of the perturbed Einstein equation ( |15[ ) gives two independent equations: 
3A£' 



(21) 



v 2 © 



(H 2 + K,)-m 

3Af 



4a 2 a 

-6H' - 12H 2 - 6/C 
3X 2 C 



2a 2 a 

(H 2 + 1C) 2 S^> -\ 
Ibcrot 2 



e' 

H(2H' - 



n 2 + K.) 



3A 2 g' 2 
8a 4 a 2 



(W 2 + AC)W 







(H 2 + K) 



3A , 



3A 
4a 2 a 



(3W' + 27i: 2 + 2/C) 



3f/2 



3A^_ 

32a 6 a 2 



(W 2 + £)W 



0' 



4a 2 c 



■{AH' + 3W 2 + 3/C) 



e 



3A 2 g' 
16a 4 a 



(H 2 +IC){2H' +H 2 



where a is 



a = 1 - ^Attf . 



(22) 



(23) 



(24) 



Eliminating Sip^ 91 " 1 from (E3) and (p3[), we can obtain one wave function for 0. However, for finding 
numerical solutions it is more convenient to directly integrate (|22] ) and (^j|) . Decomposing into Fourier 
mode and using (^) for £, we obtain the numerical solutions shown in fig. 2. Although the amplitude of 
the perturbation is somewhat enhanced near the peak of Hubble parameter, there is no growing mode 
appearing in the metric perturbation. 

Fig. 3 shows the effective density contrast defined by: 



Sp 
P 



5 s G°q 



(25) 



We notice that the small scale effective density contrast grows in Friedmann phase. In late time in 
Friedmann phase, the effect of the Gauss-Bonnet term becomes negligible and the evolution of the density 
contrast can be considered as in the ordinary Friedmann universe. Thus the growth of the density contrast 
is the consequence of the free scalar field remaining in the Friedmann universe. This result is analytically 
understood as follows. Putting A = in equations (03) and (23), we have a equation for 0: 



0" - V 2 + 6W0' + (10W + 2OW 2 )0 = 0. 
Using back ground equations (0)(||)@ with A = 0, and introducing 

a0 

u = — 

ip' 

H 



atp' 



(26) 

(27) 
(28) 



the equation can be written in a simple form[[14l: 



V 2 u- — u = 



(29) 



Regarding the perturbation as plane wave, the solution for this equation can be easily found in both 
large and small wave number limits. For large enough wave number (k 2 ^> 6" /&), the solution of the 
equation^) is 

u oc exp(ztifcry). 
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Figure 2 
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Behavior of scalar perturbation. Perturbations with different wave numbers are shown here, 
is no growing mode. 

Figure 3 
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Figure 3: Density contrast calculated using the solution shown in Fig. 2. Small scale density contrasts 
grow as y k t 1 / 3 due to the free scalar field in the Fricdmann universe. 



5 



Transforming this back into the gauge-invariant perturbative variables, for density contrast we have 

5 -l cx a cx jy 1 / 2 cx t 1 / 3 . (30) 
P 



Thus the density contrast grows like i 1 / 3 . 
3.3 Vector part 

The perturbed metric for vector perturbation is defined as: 



^ = « 2 W( _° S . Fy + Fw ) (31) 

Vector function Si and Fi satisfy S^" 1 — Fi' 1 — 0, and the combination Pi — St + F[ is unchanged under 
the gauge transformation. Vector part of the perturbed Einstein-like equation ( |l5| ) gives 2 non-trivial 
equations: 

(l-^m')(^ 2 + fC)Pi=0, (32) 

{(! - A n ?)( p i\i + P o\i)V + W - A^')(% + Pj\i) = 0. (33) 

Equation ( |32| ) shows the spatial behavior of the vector perturbation, indicating the fluctuation of curva- 
ture scale. Equation ( j33|) can be easily integrated to give 

P A3 + Pj\i cx a-\\ - ^mr 1 = (34) 
which indicates that the vector perturbation decreases as the universe expands. 

3.4 Tensor part 

The tensor part of the metric perturbation is defined as: 

^ = *0l)( J £ ), (35) 

which satisfies the constraint h*i — 0, = 0. Tensor part of ([l5]) gives one non-trivial equation for i-j 
part: 

{(1 - ^a«OftV}' + 2 ^(! ~ a^^')^' ~ (1 + ~ ^e")(V 2 - - 

Expanding ft,^ with transverse-traceless basis tensors as 

hij =h+(k,r))e +ij (k) + h x (k,r))e xi j(k). (36) 
for each polarization mode the function h(k) satisfies the equation of motion using the physical time: 

h+{3H+-)h+^-^-{l-^)h = (37) 

Numerical solutions for h(k) with different wave numbers ,in the same background as scalar and vector 
perturbations, are shown in the fig. 4. There is an instability in the de Sitter-like phase, which results 
from large positive value of | in the equation (|3^). 



(i 



Figure 4 
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Figure 4: Tensor perturbation in the same background. Note that the tensor perturbations are plotted 
with logarithmic scale. Growing mode appears in the initial de Sitter-like phase. The smaller the 
perturbation scale is, the larger the growth rate becomes. 



4 Conclusion 

We have made a perturbative analysis of the cosmological model presented by Rizos and Tamvakis |l3| , 
to investigate the stability of the model. The perturbative equations of motion are solved numerically, 
and we found that the system is unstable under tensor perturbation. This instability appears in the de 
Sitter-like stage, and kinetic-driven inflation is thought to be responsible for the growing mode of the 
perturbation. Since the form of the function £ coming from superstring effective action gives different 
asymptotic behavior of H, from this calculation alone we cannot conclude that this instability continues 
from the infinite past. But anyway, even with the modulus-geometry coupling of the form (||), there 
exists a tensor-part instability at least right before the Hubble parameter peak. Taking this model more 
seriously, we have to point out some problems. First of all, we have to check the validity of the effective 
action (||). Since the inclusion of the 1-loop effect predicts somewhat different character of the solution, 
it is certainly possible that the higher loop correction might change the nature of the cosmological model 
drastically. In this analysis we have ignored the effects of dilaton, but of course we must take these into 
acount. Viewing this model as a realistic model of our universe, we also have to find a mechanism to 
create the ordinary matter, i.e. there have to be something like reheating in the inflationary universe. 
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